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ABSTRACT 

We review the most decisive currently available measurements of the surface effective temperatures, Teff, of 
white dwarf (WD) primaries in cataclysmic variables (CVs) during accretion quiescence, and use these as a 
diagnostic for their time averaged accretion rate, (M). Using time-dependent calculations of the WD envelope, 
we investigate the sensitivity of the quiescent Tgff to long term variations in the accretion rate. We find that 
the quiescent T^ff provides one of the best available tests of predictions for the angular momentum loss and 
resultant mass transfer rates which govern the evolution of CVs. While gravitational radiation is completely 
sufficient to explain the (M) of strongly magnetic CVs at all Porb, faster angular momentum loss is required 
to explain the temperatures of dwarf nova primaries (non-magnetic systems). This provides evidence that 
a normal stellar magnetic field structure near the secondary, providing for wind launching and attachment, 
is essential for the enhanced braking mechanism to work, directly supporting the well-known stellar wind 
braking hypothesis. The contrast in (M) is most prominent for orbital periods Porb > 3 hours, above the so- 
called period gap, where (M) differs by orders of magnitude, but a modest enhancement is also present at 
shorter Porb- The averaging time which (M) reflects depends on (M) itself, being as much as 10^ years for 
low-(M) systems and as little as lO-' years for high-(M) systems. We discuss in some detail the security 
of conclusions drawn about the CV population in light of these time scales and our necessarily incomplete 
sample of systems, finding that, due to the time necessary for the quiescent T^ff to adjust, the consistency of 
measurements between different systems places significant constraints on possible long-timescale variation in 
M. Measurements for non-magnetic systems above the period gap fall below predictions from traditional stellar 
wind braking prescriptions, but above more recent predictions with somewhat weaker angular momentum loss. 
We also discuss the apparently high Teft's found in the VY Scl stars, showing that these most likely indicate 
(M) in this subclass even larger than predicted by stellar wind braking. 

Subject headings: binaries: close — novae, cataclysmic variables- stars: dwarf novae — white dwarfs 



1. INTRODUCTION 

The evolution of short-period binaries containing a Roche 
lobe filling low-mass main sequence (MS) star transferring 
matter onto a white dwarf (WD), the bulk of the cataclysmic 
variables (CVs; Warner 1995), has long been believed to be 
driven by two distinct angular momentum loss mechanisms. 
Each of these has a distinct rate at which angular momen- 
tum is lost, J, and resulting time-averaged mass transfer rate 
(M) required such that the MS star remains within the Roche 
lobe dimensions. The absolute minimum J for the binary is 
set by losses due to gravitational radiation from the orbital 
motion itself (Faulkner* 1 97 ItlPaczynski & Sienkiewicd 19811: 
iRappaport et al. 1982). Higher J is obtained due to the partial 
magnetic attachment of stellar winds similar to t hose which 
spin down isolated stars (IVerbunt & Zwaanlll981h . which ex- 
tract angular momentum from the orbit because the Roche 
lobe filling MS star is tidally locked to t h e orbit. It was 
suggested (Paczynski & Sien kiewicd 119831: IRappaport et al.l 
1 19831: iSpruit&Ritteriil983() that the evolution of CVs could 
be explained by combining these two mechanisms: at long 
orbital periods, Po^b ^ 3 hours, when the star has mass Mms 
O.25M0, and can support a typical stellar magnetosphere. 



wind losses dominate /, but when Mms S O.25M0, the star's 
magnetic structure is disrupted by the loss of the radiative core 
so that J falls to the level of gravitational radiation. 

This interrupted magnetic braking (1MB) scenario (see e.g. 
iHameurv et al.lfT988l:lKoibTT99l:lHoweU et al.l2001l) explains 
an essential feature of the observed CV population: the lack 
of systems with 2 hours < Pgy^, < 3 hours. Systems in this 
range are predicted to be detached and non-mass transfer- 
ring because the high mass loss rate to which the MS star 
was subject during the wind J phase is large enough for it 
to become bloated, such that when J decreases and it re- 
turns to the equilibrium radius for a main sequence star of 
the appropriate mass, it is well within the Roche lobe. The 
orbit must then contract further, to shorter Pmb, before con- 
tact and thus mass transfer is reestablished. Such bloating 
above the period gap does appear to be borne out by recent 
studies (Beuerm ann et al . 1998; Knigge 2006). While the re- 
lation derived from the spindown of cluster stars can, with 
questionable justification, be extrapolated to the spins appro- 
priate to the CV Porb values, what has actually taken place 
in practice is that several variations were considered early on 
(Rappaport et al. 1983) and that which best reproduced the 
observed period gap was chosen. Thus the mass transfer rate. 
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M, which 1MB predicts above the period gap is that which 
is necessary to bloat the MS star by the necessary amount to 
reproduce the Porb extent of the gap itself. 

While the 1MB scenario of CV evolution does explain the 
observed period gap, it faces a number of problems. On one 
hand, observations of single low-mass stars do not show ev- 
idence for a change in spin-down rate at t he mass bound- 
ary w here the stars become fully convective jAndronov et alj 
12003b . On the other hand, most attempts to observationally 
measure M in both J regimes h ave proven difficult. Initial 
indicators from disk properties dPattersonl [T984) were con- 
sistent with the 1MB picture, but they lack accuracy as they 
rely on a number of poorly known system properties, includ- 
ing their distances. Even if successful, M determined from 
accretion disk studies measures only the instantaneous mass 
transfer rate. Indicators for (M) from the departure from ther- 
mal equilibrium in the donor star appear quite promising, but 
suffer from the fact that reproducing the radii of isolated stars 
with these theoretical models is non-trivial (Beuermann 2006; 
iRibas 2006), and the value of M derived can be very sensitive 
to the radius predicted by the model. Townsley & Bildstenl 
(l2005h were able to use the period-specific classical nova rate 
to investigate the M-PqA relation since the nova ignition mass, 
and therefore inter-outburst time, depends on (M). However 
absolute M measurements with this method are elusive due 
to uncertainties in the mass distribution and CV population 
density. Measurements of the WD effective temperature, T^a, 
the subject of this paper, provide one of the more promising 
avenues for constraining the M in observed systems. The qui- 
escent Teff expected from a given time average accretion rate 
(M) can be calculated dire ctly (iTownsley & Bildstenl 120031 : 
iTownslev & Bildstenl llOOl herafter TB03 and TB04), and 
our discussion will largely be concerned with how represen- 
tative the (M) inferred %s is of the true (M). 

Given that the WDs in CVs are relatively hot objects, 
Tes ^ lOOOOK, their spectral energy distribution peaks in 
the ultraviolet (UV), and it is in this wavelength range that 
most of the CVWD Teff measurements were obtained. A 
handful of bright CVWDs were intensively studies with 
the Intern ational Ultrav iolet Explore r (lUE), e.g. VW Hyi 
jMateo & Szkodvl [Tosl. WZ See (ISion et a .1 fT990h . or 
AM Her ( Heise & Verbuntlll988HGansicke et al II1995I) . Tem- 
perature estimates obtained prior to the launch Hub ble Space 
Telescope (HST) were summarized by 'Sion' ('1991'), and up- 
dated for measurements obtained predominantly with the HST 
first-generation UV spectrographs by SionI ( Il999l) . A sig- 
nificant increase in the number of reliable CVWD Teff mea- 
surements has become available since the deployment of 
the Space Telescope Imaging Spectrograph (STIS) on HST 
and the launch of the Far Ultraviolet Spectroscopic Explorer 
(FUSE), see e.g. lArauio-Betancor et all ( l2005bl) . 

We begin by reviewing the relationship between (M) and 
Teff as discussed by TB03. The amount of variation expected 
in Teff due to long term M variations is evaluated from both 
quasi-static models and time-dependent envelope simulations. 
Following this, in ^ features of the system such as outburst 
properties, accretion geometry, and possible additional heat 
sources are discussed to evaluate their impact on the utility of 
Teff as and indicator of (M). ^critically reviews the available 
measurements of CV WD Teff's, with the intention of select- 
ing a well-understood set of measurements from which firm 
conclusions on CV properties can be drawn. After some dis- 



cussion of general uncertainties and biases, ^ sets out our 
conclusions, including the contrast in (M) across the gap, the 
enhanced (M) in VY Scl-type novalikes displaying low states, 
and the case for the suppression of wind braking in CVs with 
highly magnetic WDs. 

2. RELATIONSHIP OF Teff TO M 

The luminosity streaming up through the surface of the WD 
during accretion quiescence is released at all depths effec- 
tively down to the base of the accreted layer by the com- 
pression of material as it is pushed deeper into the star by 
further accretion (TB04). The infall energy is deposited at 
very shallow depths and therefore radiates away quickly 
hours) after the cessation of active accretion (see ^2. 3 1 below 
for more discussion). Due to the lengthening thermal time 
with increased depth, a large portion of this quiescent lumi- 
nosity reflects an accretion history averaged over timescales 
which can be longer than 10"* yr, depending on the character- 
istic value of M. 

As a reference point from which to begin, ^2.1 [ reviews how 
the quiescent luminosity depends on the time averaged value 
of the accretion rate, (M), and other features of the WD. This 
provides a relation between Teff and (M) which has no "free" 
parameters, but does have important but manageable uncer- 
tainties related to unknown properties of the WD. Follow- 
ing this, in ^2.21 we present a basic discussion of how the 
response of the envelope to accretion can be understood in 
terms of the run of local thermal time with depth. This pro- 
vides justification for the general assertions made above about 
the timescales on which the surface luminosity can vary. Us- 
ing estimates of the run of thermal timescale in the outermost 
portion of the envelope, in 32.31 we derive the timescales on 
which heat deposited near the surface is radiated away. This 
demonstrates how measurement of the luminosity escaping 
due to compression is possible between transient accretion 
events such as dwarf novae. Such simple estimates are in- 
sufficient for cha racterizing long-term variations, so that we 
proceed, in 32.41 to discuss and then characterize, with sim- 
ulations, how long timescale variations in M affect the quies- 
cent Teff. 

2.1. Long Tenn Average Properties 

|TB04| presented a detailed discussion of the impact of ac- 
cretion at rates low enough that hydrogen burns in unstable 
outbursts. Classical Nova eruptions, on the thermal structure 
of a WD. For this paper we are primarily interested in the pre- 
dictions for the quiescent surface luminosity, Lq(M, {M),Tc), 
where M is the WD mass, (M) is the time-average accre- 
tion rate, and Tc is the WD core temperature. Through 
a simple balance of compression of material and radiative 
heat transport, the quiescent surface luminosity is given by 
Lq ~ {M)Th/ iiirip = AirR^aT^ff where 7}, is the temperature at 
the base of the radiative layers, typically Tt, ~ T^ R is the 
WD radius, p is the mean molecular weight of the accreted 
material, nip is the mass of the proton and cr is the Stefan- 
Boltzmann constant. TB04 found that Tc approaches an equi- 
librium value when the WD is subject to accretion at con- 
stant (M) for timescales similar to the WD core thermal time, 
^ 10^ years. This Te eq is not strongly sensitive to M and 
increases weakly with (M), being generally 0.5 to 1 xlO^ 
K for CV WDs, which have (M) = 10-'' to IO-'^Mq yr"'. 
This equilibrium has recently been demonstrated in Simula- 
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tions which follow the WD evolution through many nova out- 
bursts, finding slightly but not significantly higher Tc.eq for the 
higher (M) 's in this range, and similar evolutionary timescales 
(jEpelstain et al. 2007). Lq is insensitive to for T^. < eq and 
the (M) relevant above the period gap, where (M) changes on 
timescales less than the core thermal time (TB03). Thus we 
do not need to know the details of 7^. in order to derive (M) 
from Teff , though we use eq as a convenient reference value. 

Including a small amount of nuclear heating, using the 
methods of TB04 we find that (at T^^eq) the average Lq dur- 
ing the classical nova (CN) cycle is 



Lq=6x 10" 



10-»'Mq yr- 



,O.9M0y 

Using an approximate power-law relation for R oc M"'^ near 
M = 1 .OMq, this gives 



7;i=f= 1.7 X 10* K 



(M) 



IO-i'^Mq yr- 



1/4 



M 



Q.9M, 







(2) 



As emphasized in TB03, the (M) inferred from a system's 
quiescent Tgff is strongly dependent on the assumed M, almost 
entirely due to the use of R to infer Lq. While (M) cx for 
agiven Leff, (m) = {M)/4t:R^ cxM"°^ Just a25% uncertainty 
in mass, allowing M = 0.6 to I.OMq, leads to nearly a factor 
often uncertainty in (M). This seems to be the strongest limit 
on the utility of measurements of Lgff, and is difficult to avoid 
without either independent mass measurements or distances. 
Only a handful of systems have independent mass measure- 
ments (see Table [U, and all but one are DN below the period 
gap, and the one above is a VY Scl star with a high (M). It 
appears infeasible to make conclusions based on a subsample 
that only includes objects with mass measurements. However, 
as we will see below, order of magnitude contrasts are easily 
discernible when the full sample is considered and are impor- 
tant for discriminating between angular momentum loss laws. 
Additionally, there are now enough sound measurements that 
further progress can be made with some assumptions about 
the mass distribution of the population. 

The next most important source of uncertainty, which is 
nearly impossible to eliminate for most systems, is due to the 
increase in Lq as the mass of the accumulated layer, M^^^, in- 
creases between classical nova outbursts. Scatter in the ob- 
served Leff values at the level of ±5% is expected due to sys- 
tems having different Mace's at the current epoch. This degree 
of scatter is drawn from allowing a range O.OSMign < Mace < 
0.95Mign, and, though the full allowed range gives a larger 
variation, this represents well the size of layer an observed 
WD is likely have. (See TB04 for details of how Lgft varies 
with Mace-) We set the mass fraction Xiy^^ = 0.001 in the ac- 
creted material throughout this work, because the difference 
between this and similar predictions for Xi-^e = 0.005 is less 
than the uncertainty due to the unknown Mace (TB04). 

2.2. Envelope Response Timescales 

One of the attractions of utilizing Leff as an indicator of (M) 
is that Leff contains averaged information about the time his- 
tory of M rather than its instantaneous value. This is of great 
benefit for comparing (M) with predictions based on long- 
term drivers of the binary evolution contributing to orbital an- 
gular momentum loss. Before proceeding with simulations of 
the response of the envelope to (M) variations, it is useful to 



first discuss the relevant timescales within the envelope and 
how they relate to the surface luminosity. We do this in order 
to understand how much a perturbation on some prescribed 
timescale is likely to affect the surface luminosity. 

The heat equation in the outer layers of the WD is given by 
(TB04) 

dT dF T 
^''li7 = 7r + ^^""-(^-^^d), (3) 
at ay y 

Where y = J* pdr ^P/g, the column depth in from the surface 
(r = R), forms a radial coordinate, P is the pressure at radius 
r, g is the surface gravity, F = {4aT^ /3K)dT /dy is the local 
area-specific energy flux with a being the Stefan-Boltzmann 
constant, k the opacity, cp the specific heat at constant pres- 
sure, V = dlnT /dlnP, Vad = 9 In r/91nP at constant entropy, 
and »j = M/4ttR^. The temperature profile of the envelope at 
a given depth, y, can change on the thermal time for that layer, 
which we define by "one-zone" differencing the left hand side 
and the first right hand term in equation (O, dropping the ac- 
cretion source term, to obtain. 



Tth ; 



ycpT 3Ky cpT 
F ~ AaT'^ 



(4) 



In a static envelope state, which is a good approximation for 
a steady M (TB04), dT /dt is small or zero, so that dF/dy, 
the contribution to the surface flux from each layer, is set by 
m and local properties of the layer. Thus we expect that vari- 
ations in M will appear as variations in dF/dy\y after a time 
tthiy), when that layer can respond. For order unity variations 
in M on a given timescale fvar, the contribution to Lq from lay- 
ers which have tth > fvai will not be affected by the variation, 
while the contribution from layers with fth < fvar will change 
by unity along with M. Thus L(surf)-L(rth = fvar), the contri- 
bution to the surface luminosity from the layers outside that 
with fth = fvar, evaluated with the average envelope state, pro- 
vides a reasonable indication of the possible variation in Lq 
which can be expected due to unity-level variations in M on a 
timescale of fvar- This is very similar to how the characteristic 
cooling time of a dwarf nova outburst is set by t he thermal 
time at the bottom of the freshly accreted material dPiro et al] 
12005ft . 

This analysis no longer applies in the deepest layers where 
the thermal transport becomes mediated by electron conduc- 
tivity. This is because the one-zoning used to define r,h no 
longer holds. Once fvar exceeds the thermal time of the whole 
radiative layer, we are left with a problem which is more 
similar to the classic WD cooling problem: an insulating 
layer which acts as the thermal regulator for the underlying 
heat reservoir formed by layers which are much more ther- 
mally well-coupled by electron conduction. In this case the 
timescale for further changing Lq is expected to approach the 
CN inter-outburst time, that required to build up the maximum 
degenerate region. 

2.3. Reaching the quiescent Leff 

With the above estimate for thermal response time with 
depth, it is useful at this point to discuss specifically the cool- 
ing of the thin outer layer that is heated by the infall energy 
of the accreted matter This is essential to understanding why 
the quiescent Lgff is a good indicator of the energy being lib- 
erated by compression in the deeper layers of the star When 
material reaches the surface of a WD via the disk, half of the 
gravitational infall energy (GM/R per unit mass, where G is 
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Newton's constant) has been radiated in the disk and the rest is 
possessed as kinetic energy. This kinetic energy is deposited 
immediately at the WD surface as the material is stopped and 
spre ad by interaction with the surface layers (Piro & Bildsten 
l2004l) . This process can heat the WD surface to quite 
high temperatures given by r^urf = [(GMM /2R) /AnR^a] '/^ « 
10^ KiM/lQ-^MQ yr-i)i/4(M/O.9M0)i ''. 

However, this heat does not penetrate the WD due to the 
opposing thermal gradient in the underlying radiative atmo- 
sphere which has a profile characterized by Fg = Lg/AirR^ 
satisfying 4-aT'^/3Ky = Fq. During accretion, the infall heat- 
ing can penetrate to approximately where the local temper- 
ature is the same as Tsurf, or a column depth of jsurf — 
AaT^^^^/linFa^. Thus, if n is assumed to be dominated by 
electron scattering at this depth for simplicity, the mass 
of the layer heated by infall is Mjurf = 47r/?^3'surf — 5 x 
10-''Mo(7;urf/105K)4(reft;q/14kK)-4. The thermal time for 
this layer to cool is between ycpT /F^wf and ycpT /Fq which 
are 10^ and 10^ seconds for Teff^q = 14 kK. The latter is fairly 
consistent with the 2.8 days fo und for the cooUng of V W Hy i 
following a normal outburst (Gans icke & Beuermannll996al) . 
This timescale should be typical for cooling after a normal 
dwarf nova outburst, but is much shorter than the cooling time 
after a superoutburst, in which an order of magnitude more 
material than this is deposited. In that case, in contrast, heat 
is released on the cooling time f or the a dded material due to 
its relatively rapid compression (iPiro et al. 2005). Note that 
in reality much of the ma terial added to th e WD cools as it 
spreads over the surface dPiro & BildstenI 12004) so that the 
above analysis only applies to the heated region near the equa- 
tor while the rest of the star remains near T^ff q. 

2.4. Numerical Simulations of Long Term Accretion Rate 
Variation 

The analysis in ^2.21 indicates that a given variation in M 
leads to a larger variation in the quiescent surface flux, Lq, 
if it occurs on a longer timescale. That is, the excursion 
from the average quiescent luminosity SL^ = max |Lq- (Lq) | « 
(i'q,max-iq,min)/2 incrcascs with the timescale of M variation, 
fvar, for a fixed amplitude of M variation. Thus very brief 
variations in M, such as dwarf novae outbursts, have very lit- 
tle impact on Lq despite their large magnitude in M, once the 
shallow transient has passed. The luminosity excursion, (5Lq, 
of course also depends on the magnitude of the M variation. 
This dependence is expected to be fairly simple, roughly lin- 
ear, and therefore we focus here on the dependence on fyai, 
which relates to the thermal response time structure of the en- 
velope. 

In order to test how Lq, and therefore the quiescent Left, will 
vary in response to M variations, we have performed simu- 
lations of an accreting layer on a WD in the plane-parallel 
approximation. This approximation is fairly good through the 
accreted layer and has little impact on the time-variation prop- 
erties being studied here. A similar appro ximation, extend ing 
only to more shallow depths, was used bv lPiro et al.l (|2005|) in 
studying the decline of T^ff after a dwarf nova superoutburst. 
Our numerical treatment is described in the appendix. 

As simple experiments, we have applied two forms 
of long-term M variation: sinusoidal, with M = (M)(l - 
0.8sin[7rf/fi/2]), and square wave, with M = (M)(l - 
0.8 sign[sin(7rf/Pi/2)]), for a wide range of half-period of M 



variation, Pi/2- Short term variations in M such as accretion 
disk cycles (dwarf nova outbursts) are not treated explicitly in 
order to keep timesteps large. Figure [T] shows time series of 
Lq and Leff for variations in M on two timescales, P1/2 = lO'' 
and 4 x 10"^ yr This example uses M = 0.9Mq and (M) = 5 x 
1O~"M0 yr"', representative of DN systems below the period 
gap, and has Leff^q = 14 kK and Mig„ = 2.3 x 10"'*Mo. Time 
zero corresponds to Ma^c = O.SMign. The model is initialized 
by starting from the static solution described in TB04 with 
Mace = 0.25Mign and applying M = (M) until Mace = 0.5Mign to 
allow the model to settle on the new numerical grid. In Fig- 
ure[ri a modest increase in the total variation is observed with 
the square wave M applied instead of the sinusoid. There is 
no appreciable phase difference between the applied M vari- 
ation and that the response. As expected, for the same (M) 
and magnitude of (M) variation, a longer timescale leads to 
a larger variation in Left. Consecutive peaks are slightly in- 
creasing due to the increasing Mace during the buildup to clas- 
sical nova, the longer displayed example making it to nearly 

Mace = 0.59Mig„. 

There is a tremendous contrast between observable 
timescales of tens of years, and the time between classical 
nova outbursts which can be as much as 10^ years depending 
on (M). In order to probe this variety, we have performed 
the same 5-cycle simulations shown in Figure[T]for variability 
timescales up to about 2.5% of the classical nova accumu- 
lation time, making the total simulation time about 25% of 
the accumulation time. We characterize the variation in Lq 
by its excursion (5Lq k, (Lq.max -i'q,min)/2, where the maxi- 
mum and minimum are evaluated for each cycle and the dif- 
ference is then averaged over all cycles. This is then divided 
by the mean, (Lq), to obtain a fractional variation. Figure |2] 
shows the excursion found from the numerical simulations us- 
ing square (solid lines) and sine (dashed lines) wave M. Three 
cases are shown, M = O.9M0 with (M) = 2 x lO-'^M© yr"' and 
M = 0.9 and O.6M0 with (M) = 5 x 10"' 'Mq yr"' . These have 
(iq) (7;ff) of 120, 2.8, and 2.6 x IO'^Lq (35, 14, and 12 kK) 
respectively, Mig,, of 2.2, 23 and 39xlO~^M0, time between 

classical novae of 1.1 x 10"^, 4.6 x 10^ and 7.8 x lO*" yr, and 
L, of9.6,5.7and5.7x 10'' K. 

As discussed in ^2.21 the magnitude of the variation or ex- 
cursion in Lq can be estimated by evaluating the compres- 
sional heat release at depths with thermal time less than fyai. 
This can be expressed as L(surf)-L(fth = fvar), where the qui- 
escent surface luminosity, L(surf) and the luminosity at the 
depth where the thermal time matches the variation time of 
the applied M variation, L(fth = fvar) are evaluated from the 
run of L and fth in a static accreting envelope at M = (M), 
L(surf) = (Lq) and Mace = 0.5Mign. These curves are shown 
for comparison in Figure |2] one for each case shown from the 
time-dependent simulations. 

We find that L(surf)-L(fth = fyai) provides a good estimate 
of the variability until electron conduction becomes impor- 
tant, fth ^10"^ yr. Deeper than this, the fth we have defined 
underestimates the thermal time, so that L(fth = fyar) is being 
evaluated at too deep a layer and L(surf)-L(fth = fyar) is there- 
fore being overestimated. The variability seen in Lq in the 
simulations does not increase for longer timescales, due to 
the increased heat capacity from the thermally connected lay- 
ers where electron conduction provides heat transport that is 
efficient relative to the overlying radiative layer We find that 
the averaging time for a given system depends strongly on its 
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Fig. 1 . — Time evolution of quiescent surface flux. Lq = AirR^aT^ff, due to compressional energy release in the envelope for a WD with M = 0.9Mq and 
(M) = 5 X 1O""M0 yr"'. Two variations are applied, a sine (dashed lines) and a square wave (solid lines) both with amplitude of 0.8(M). Two timescales of 
1/2 = (Isft) and 4 X lO"* yr (right). For the same magnitude variation in (M), a longer timescale of variation leads to a slightly 
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Fig. 2. — Fractional excursion in quiescent surface flux, SLq/{Lq) = 
(i-q.max ~ i'q.min )/2 {Lq ) , as a function of timescale of the applied M variabil- 
ity, ?var. The response of an envelope simulation to a 0.8(M) variation in M 
is shown for a square (solid) and sinusoidal (dashed) variation. One mass is 
shown for (M) = 2 X 10""' yr"', with M = 0.9Mq, (far upper) and two 
for (M) = 5 X 10-" Mq yr-', with M = 0.9 (middle) and 0.6Mq (lower). 
These cases have (Lq) {(T^fi)) of 120, 2.8, and 2.6 X IO'^Lq (35, 14, and 12 
kK) respectively. Also shown for comparison (dot-dash lines) is the run of 
[L(surf)— L(f,i, = ?var)]/ {Lq) in a static approximation of the average envelope 
state. This represents at what depth, in terms of timescale, energy is being 
released by compression, and provides a good estimate of the variation ex- 
pected in response to order unity M variation except at the deepest parts of 
the envelope where electron conduction becomes important. 



(M). This is due to the fact that while the luminosity varies by 
two orders of magnitude (approximately oc (M)) for CVs, the 
thermal content of the radiative layers of the envelope changes 
little with (M) . Note that we have only explored order unity M 
variations here for simplicity of demonstration, smaller varia- 
tions will lead to correspondingly smaller responses. 

The magnitude of the variation in quiescent surface lumi- 
nosity, Lq, depends on the timescale of the proposed varia- 
tion in M. Since the the (M) inferred from Teff is approxi- 
mately proportional to the corresponding Lq, fractional varia- 



tions induced in Lq by a time-variable M correspond directly 
to fractional uncertainty in an inferred (M). Since Figure |2] 
shows the response to an order-unity variation in M, it also 
approximately quantifies the fractional uncertainty for a given 
timescale of variation. Assuming that M has only short or 
moderate timescale variations, even up to lO-' yr, but is consis- 
tent on longer timescales, Leff is a good indicator of the time- 
averaged M, with minimal uncertainty (20%) at low M and 
moderate (50%) at highM. However, due to the long but finite 
thermal time of the envelope, under the assumption that M 
varies over very long timescales (> 10^ yr for high-M systems 
or > 10^ yr for low-M systems), Leff is not a reliable indicator 
of (M), having a large uncertainty. In this case T^ff depends 
more strongly on the recent accretion history than the over- 
all average. Thus the character of the mass-transfer law being 
tested must be considered when Leff measurements are used to 
infer (M) values. The utility of this relationship between un- 
certainty and timescale and its application to candidate mass- 
transfer scenarios is discussed in more detail in ^ There, 
multiple Leff measurements from similar systems are used to 
obtain an independent constraint on the time-variability of M, 
removing the a priori model dependence for some systems. 

3. SYSTEM FEATURES AFFECTING THE UTILITY OF Teff 

In addition to long-timescale variations in the mass trans- 
fer rate, several features of the systems in which the quiescent 
Leff measurements can be made affect our ability to determine 
Lq from observations. While we believe that these issues can 
be avoided by careful selection and interpretation of observa- 
tions, it is worth summariz ing the relevant issues to justify this 
assertion. iPiro et all (l2005i) have performed an excellent anal- 
ysis of how a dwarf nova cools after outburst, finding that the 
cooling time depends on the amount of matter accreted in the 
outburst. We therefore forego a detailed discussion of those 
systems, only noting that care is typically taken that Left is 
measured as far in quiescence as is feasible. The other impor- 
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Fig. 3. — Decline of T^ff for the 3-year low state of the novalike TT 
Ari during 1982-4. UV spectral observations taken at approximately 1 and 
2 years after the decline into the low state are both consistent with T^if = 
39000 ± 2000 K ( Gansicke et al. 1999). The three example cases shown have 
the same (M) = 4.1 X IO'^Mq yr"' but different high state durations: 30 yr 
(solid), 10 yr (dashed), and 90 yr (dot-dashed). 

tant kinds of systems in which Teff can be measured appropri- 
ately are novahke variables during extended quiescence inter- 
vals and Magnetic CVs, again in quiescence, where material 
impacts on a small portion of the WD surface, and the emis- 
sion from the rest of the star, reflecting Lq, can be separated. 
After discussing these two cases we provide a direct estimate 
of the timescale on which a CVWD would cool to Teff q after 
a thermonuclear runaway, and close with a brief discussion 
of additional energy sources related to accretion which could 
cause the Tgff observed to differ from our calculated Lq. 

3.1. Cooling From High State 

In the VY Scl stars discussed in section l4~2l the high M ^ 
\Q~^Mq yr"' will, on rare occasions, turn off for timescales 
of up to a few years. This provides a windo w in which to 
measure the T^f^. As discussed above in section l272l the layer 
which is heated by infall cools in ~100 seconds at this high 
M and resulting high Teff q. Since the outer layers, which can 
change their thermal state during quiescence, contribute a rel- 
atively small fraction of the overall Lq, the cooling during the 
quiescence is modest. In Figure |3] we have calculated an ex- 
ample evolution for the ~3-ye ar quiescence of TT A ri during 
which it was observed twice (iGansicke et alJ [19991) . assum- 
ing various high state durations between regularly repeated 
quiescent intervals. With this type of analysis we are able to 
solve directly for the actual (M) = 4.1 ±0.8 x lO'^M© yr"' 
forM = O.9M0, accounting for the cooling of the outer layers 
during quiescence. It is observed from Figure |3] that the Teff 
flattens out as the thermal time of the cooling layer becomes 
longer, in such a way that the two epochs of measurements 
(shown) will be the same within the observational error By 
varying the high state interval at fixed (M), we see that, even 
for durations as short as 10 years, the quiescent flux is a good 
indicator of (M) instead of the high-state M. The Lq scales 
linearly with (M) as expected. 

3.2. Polar Accretion 

The geometry of the accretion on the WD surface differs 
significantly for a strongly magnetic WD {B ^ 10^-10^ G), in 
which case mass is deposited at the magnetic poles. How- 
ever, due to the depth at which this material must spread over 
the star and the depths at which energy is liberated by com- 
pression, the effect on the T^ff measured away from the polar 
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Fig. 4. — Luminosity as a function of depth as measured by P as a fraction 
of Lsurf. The three cases shown are for (M) = 2 X 10"' (lower) and 5 X 
1O""M0 yr"' (upper) at M = 0.6Mq (solid) and only the lower (M) at M = 
0.9Mq (dashed). These have (reff) of 140, 2.6 and 2.8 X 10"^ Lq (32, 12, 
and 14 kK). The bulk of the compressional luminosity is released below the 
depth at which spreading occurs. 

regions is modest. The magnetic field can only constrain the 
accreted material to the polar regions down to a limited pres- 
sure depth, Pcrit, below which lateral pressure gradients are 
strong enough to force the material to spread over the surface. 
Heat liberated by compression up to Pent is localized to the po- 
lar regions, while that from compression at higher pressures 
is spread over the whole surface. 

The magnetic field can keep the accreted material con- 
strained to the WD poles up to a depth where /J^rit — 21 /h 
(iHameury et al.lfT983h where (3 = P/A-iiB^ is the ratio of gas to 
magnetic field pressure, I is the size of the polar cap, and h 
is its thickness. Using h = P/ pg, t = 10^ cm and a degenerate 
equation of state we find that P^nt ^ lO'^ erg cm"^ (gshB^)^/^ , 
where a subscript s indicate values divided by 10' cgs units, 
e.g. gs = g/lQ^ cm s^. Note that a nondegenerate calcula- 
tion, which is more appropriate since degeneracy sets in at 
P ~ 10'^ erg cm"^ for Tc= 10^ K, would give even lower Pent- 
We forgo such a calculation here since the difference does not 
impact our conclusions. 

As discussed above, energy is liberated by compression at 
all depths, with much coming from deep in the accreted layer 
Figure ID shows how L/Lq varies with depth into the star as 
measured by P. In a magnetic CV WD, we assume that ac- 
creted matter spreads over the WD surface at approximately 
Pcrit, so that only L|p=p„j, is released over the whole surface 
of the star For typical fields of 10^ G, this amounts to better 
than 80% of the Lq expected in the non-magnetic case. Thus 
magnetic CV Leff's warrant a small upward correction from 
the simple relation given by eq. ^ to obtain the actual (M) . 

3.3. Cooling After Thermonuclear Runaway 

After shutoff of the nuclear burning in a classical nova, due 
to heat leftover from the outburst it will take some period of 
time for Lett to return to that characteristic of (M). Although 
this timescale could vary widely depending on how much of 
the burning envelope is actually left behind on the WD, it ap- 
pears that in most cases the ejected matter is similar to or 
larger than the total accreted, indicating that this layer is fairly 
thin (TB04). The enhancement of the temperature below the 
burning layer is expected to be modest, being limited to just 
above 10^ K by the onset of the instability. A useful estimate 
for the cooling time is then to take the cooling time of the 
outer radiative layer Using Lq, the endpoint luminosity of the 
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decline, to estimate this cooling time gives a good upper limit 
on the decline from the hot outburst state. 

The mass coordinate, measured in from the surface, where 
the free-free opacity and that due to conduction are equal is 
AM ~ 4 X IQ~^Mq{M /Q.9Mq)~'^-^ for an interior temperature 
of lO' K. This leads to a thermal time of the radiative layer of 

r,!, ^ 7 X 10^ yr V7 ' , (5) 

where we have used the approximate dependence R oc M"' -^. 
Note that this is similar to the timescale at which the re- 
sponse curves from above the variation analysis (see Figure 
121) flatten o ff. This is rea sonably consistent with the timescale 
found by P rialnild (Il986h from nova simulations at high mass, 
M = 1.25Mo. This is generally much shorter than the inter- 
outburst time and therefore we expect few objects to show 
such an enhancement. If the heated region extends deep 
enough for some degenerate material to have a longer cooling 
time, which may be the case at low (M), it would introduce an 
enhanced lower boundary temperature for the fresh envelope 
for a longer time. However, this will not lead to significant en- 
hancement of Lq because Lq is much larger than the cooling 
luminosity set only by the modestly increased temperature at 
the base of the layer. 

3.4. Possible Additional Heat Sources 

If there is a heat source depositing excess energy at the WD 
surface during quiescence, this can lead to misinterpretation 
of the quiescent Teff as representative of Lq. The energy re- 
lease due to infall is much larger than that due to compression 
in the WD, so that a small amount of steady accretion dur- 
ing quiescence in non-magnetic systems could compete with 
the flux from below in heating the WD surface. Most of the 
energy from quiescent accretion is expected to appear in the 
X-rays and typical X-ray luminosities from dwarf novae are 
1030 _ 1031 erg s-i (IVerbuntet al. l fT997h requirmg a quiescent 
M ^ 10"'^Mo yr~'. While this is clearly less than the lu- 
minosity of the brighter systems we will consider, it is only 
slightly less than the lowest Lq ~ 10^' erg s"' {T^ff = 14 kK, 
M = Q.9Mq) that we are considering here. We will ignore 
this contribution in the first approximation for two reasons. 
First, it is not expected that Luv and Lx should be similar, if 
anything Luv should be much lower, if there is a steady state 
corona above the WD surface where the infall energy is re- 
leased. Any heating which modifies Left must be deposited 
at the photosphere, and no higher. Second, while there is a 
wide variety of X-ray luminosities for dim DN which has no 
clear cause, we will show that the Leff for these systems show 
relatively little scatter At worst our dimmest systems provide 
upper limits on the (M) if some portion of the flux is attributed 
to accretion in quiescence. 

4. MEASUREMENTS OF CV PRIMARY Trff 

The most commonly used method to determine WD effec- 
tive temperatures in CVs is to fit synthetic spectra to optical 
and/or ultraviolet (UV) observations. While for single white 
dwarfs accurate temperatures and good estimates of the sur- 
face gravity (and hence, adopting a mass-radius relationship, 
the WD mass) are routinely derived from spectral fits to the 
Balmer lines alone, a number of caveats have to be consid- 
ered when applying such spectral fits to the observations of 
CVWDs. One main disadvantage encountered in CVs is that 
their optical light is a mixture of emission from the accretion 



disk or stream, hot spots on the disk edge or WD surface, from 
the donor star, and finally from the accreting WD. In the ma- 
jority of CVs, the WD contributes only a small fraction to the 
emission at optical wavelengths, which foibles any attempt 
to determine its properties from ground-based observations. 
Even in the cases where the WD is a significant source of opti- 
cal flux, implying low prevailing accretion rates, temperature 
measurements are subject to a major ambiguity: the strength 
of the Balmer lines reaches a maximum around 15 000 K, with 
the exact value being a function of surface gravity, and nearly 
equally good fits can be achieved on the "hot" and the "cold" 
side. 

Because the WDs in CVs are moderately hot, 10000- 
50000 K, their spectral energy distribution peaks in the UV, 
and therefore the most reliable information on CVWDs is ob- 
tained from space-based telescopes such as the lUE, HST, or 
FUSE. Besides a larger, or dominant, WD contribution to the 
total UV flux from the CV compared to the optical, the degen- 
eracy between hot and cold Balmer line fits is broken by the 
opacity of quasimolecular Hj and H2, which causes broad ab- 
sorption lines near 1400 A and 1600 A fo r temperatures below 
- 1 8 000 K and - 1 3 000 K, respectively dKoester et al.ll985l) . 
Examples of the degree of uncertainty in temperature s based 
on optical d ata are GW Lib and LL And, where Szkod vet al.l 
dlOOOb) and ISzkodv eTaTI (l2000ah estimated Leff ~ llOOOK 
from optical data obtained during quiescence, whereas the 
temperatures determ ined from HST/STIS spectroscopy ar e 
14000- 15000K (S zkodv et al.ll2002at iHoweU et al.ll2002h . 
While UV spectroscopy can provide fairly accurate WD ef- 
fective temperatures, its diagnostic potential for determining 
the surface gravity, \ogg, is very limited, and so is, therefore, 
the ability to derive the white dwarf mass from spectral fits. 

A remaining issue in modelling the UV data of accreting 
white dwarfs is that a large fraction display a second contin- 
uum flux component, that typically contributes 10-30%. A 
number of suggestions for the nature of this component have 
been made , such as optically thick hot accretion belts on the 
WD (e.g. lLongetal.l[T993[lGansicke& Beuermannlll996bl: 
iHuang et al.ll 199611 optically thick emission from the hot spot 
or optically thin emission from a chromosphere on the accre- 
tion disk (e.g. Gansicke et al. 2005). While the exact nature 
of this additional component is not clear, and may well dif- 
fer among the objects, its impact on the effective temperature 
determination appears to be only modest. 

An alternative method that yields WD effective temper- 
atures, as well as potentially radii and masses as well, is 
the modelling of the W D ingress/egress in mu lti-color light 
curves of eclipsing CVs (! Wood & Hornelfr990 h. The applica- 
tion of that method has been limited until recently to a hand- 
ful of bright CVs, but since the fast triple-beam CCD camera 
ULTRACAM (DhiHon ej al. 2007) became available on 4 m 
and 8 m telescopes, a number of detailed CV WD studies have 
been carried out (e.g. iLittlefair et al.ll20b6bl) . 

The literature is humming with values of CVWD temper- 
atures and some care has to be taken to differentiate be- 
tween measurements and estimates. For some purposes, it 
may be desirable to maximize t he number of available Teff 
values, such as prov ided e.g. by IWinter & SionI (l2003l) and 
lUrban & SionI (l2006h . In the context of providing a stringent 
test of the theory angular momentum loss in CVs, we focus 
here on the most accurate Lgff measurements. Table[T] lists 
CVWD Left values that we consider reliable on the basis that 
the WD has unambiguously been detected either spectroscop- 
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ically, or, in the case of eclipsing systems, through its eclipse 
ingress/egress. Below in Sect. 14. 11 - 14.31 we discuss particular 
issues that relate to the effective temperature measurements in 
the three major CV subclasses, and summarize in Sect. l4.4l 

4.1. Dwarf novae 

Dwarf novae are a subset of non-magnetic CVs with low 
mass transfer rates. The accretion disks in these systems 
are thermally unstable, and undergo outbursts lasting a few 
days to a few months with recurrence times of a few weeks 
to tens of years. Despite their relatively low mass trans- 
fer rates, only a relatively small number of dwarf novae 
reveal their accreting white dwarfs at opti cal wavelengths, 
e.g. W ZSge (Greenstein 1957) or VWHyi ( Mateo & Szkodyl 
[1984; Smith et al. 2006). Moving to the UV, about a third of 
all short-period dwarf novae are dominated by emission from 
the white dwarf. It is currently unclear why dwarf novae with 
nearly identical orbital periods and outburst frequencies differ 
radically in the characteristics of their UV spectra, such as e.g. 
VW Hyi and WX Hyi, wit h the first one bein g one of the best- 
studied CVWDs (iMat eo & Szkody 1984; Sion et all [1995b|; 
Gansicke & Beuerma nn 1996a; Lo ng et alJ fr996; Smit h et aP 
20061) and WX Hvi. where no convinc ing spectros copic signa- 
ture from the WD has been detected dLong et al.l l2005). The 
fraction of dwarf novae above the period gap where the WD is 
clearly discernible in the UV is much smaller than at short or- 
bital periods, as the accretion disks in these systems are larger 
and can sustain higher accretion rates while still remaining in 
quiescence. 

During dwarf nova outbursts, the accretion rate onto the 
white dwarf increases by several orders of magnitude with re- 
spect to the quiescent value, causing a short-terrn heating ef- 
fect (lSionlll995l;lGans icke & Beue rrnanrill l996bl; IC heng et alj 
I200d iGodon et all 12004; Piro et al.||2005t). Therefore, when 
determining the secular WD effective temperature, care has 
to be taken to observe dwarf novae as long as possible after 
an outburst. In systems with very short outburst recurrence 
times, it may b e that the WD never cool s to its secular tem- 
perature (G ansicke & Beuermannlll996bl) . (See also the dis- 
cussion in ^2.31 ) 

A final caveat relates to the spectral modelling of high- 
inclination dwarf novae, where the line-of-sight passes 
through absorbing material located above the accretion disk, 
termed an accretion veil, affecting the effective temperature 
determination. In the case of OYCar, a fit ignoring the 
veiling component yields T^s w 15 000 K (Horne et al. 199^, 
where as T^s ~ 17000K when taking the ab sorption into ac- 
count dHorne et al.ll 1994 ICheng et alJlIOOOl). While th e effect 
is noticeable at incUnations > 70° ("e.g. iLongetalJ (l2006l) . 
it is most problematic at higher inclinations where the WD 
eclipse offers an al ternative/independen t possibility of a T^s 
measurement (e.g. iLittlefair et al.ll2006ah . 

4.2. Novalikes and the case ofVYScI stars 

Novalike variables are non-magnetic CVs, mainly with pe- 
riods > 3 h, with high mass transfer rates in which the accre- 
tion disk is in a stable hot and optically thick state. Conse- 
quently, the flux of novalike variables is entirely dominated 
by the disk at optical and ultraviolet wavelengths. The only 
possibility to learn about the properties of WDs in novalike 
variables occurs if the mass transfer decreases, or turns off 
completely, so that the WD becomes visible. Novalike vari- 
ables who show such low states are called VY Scl stars, after 
the prototypical system, and are found predominantly in the 



orbital period range 3 - 4 h. The physical cause of the oc- 
currence of low states is not fully understo od, and may be 
related to star spots on the secondary star ( Livio & Pringld 
1994; Hessman et al. 2000) or irradiation driven mass trans- 
fer cycles ( Wu et al. 1995). 

Due to the rare and unpredictable nature of low states, only 
three novalike variab les have been studied at a sufficient leve l 
of detail: TT Ari (Shafter et al."1985'; 'Gansicke et al."1999b, 
DWUMa (Kniggeetal. 2000; Araujo-Betancor et al. 200% 
and MVLyr (iHoardet al.ll2004<) . In aU three systems, hot 
rwd>40000K are found. 

4.3. Polars 

Polars, or AMHerculis stars, contain strongly magnetic 
WDs. The rotation of the white dwarf is synchronized with 
the orbital period, and the formation of an accretion disk is 
suppressed. Accretion occurs via an accretion stream that 
feeds matter onto the magnetic pole of the WD. Polars enter 
low states with little or no accretion, and during these episodes 
the systems appear practically as a detached WD plus main 
sequence star. However, while the white dwarf is fully ex- 
posed during low states, its magnetic field complicates accu- 
rate Teff determinations, as the Balmer lines are subject to Zee- 
man splitting, and no accurate theory for the line profiles of 
the Zeeman components exists so far (lJordanlll992j ). Zeeman 
splitting is much weaker for the Lyman lines, and for fields 
< 30 MG the effect of the magnetic field on temperatures ob- 
tained from UV observations around Lya is relatively small. 
The polars in Table[T]have all fields < 30 MG at the primary 
accretion pole, with the exception of QS Tel and VI 043 Cen 
(both B ~ 56 MG). 

For stronger fields, even the Lyman lines become useless 
for temperature estimates, and, worse, the spectral models 
fail to even reproduce the UV/optical spectral e nergy distribu- 
tion, a n effect known from single white dwarfs (ISchmidt et alJ 
Il986l) . Consequently, the WD temperatures of high field po- 
lars are only very approxima^tively known, 1 5 000 -25 000 K 
for AR UMa dGansicke et al."2 001j) and 1 7 000-23 000 K for 
RXJ1554.2H-2721 ([Gansicke et al. i2004bl) . 

The highly asymmetric accretion geometry in polars results 
in heati ng part of the WD atmosphere around the magnetic 
pole(s) (iGansicke et al.lll995l) . As the magnetic axis of the 
WD is usually not aligned with its spin axis, the heated pole 
cap acts as a light house, causing a significant variation of 
the UV flux as a function of WD spin/orbital phase. The 
pole caps are observed also during low states, and it is not 
clear if this is due to de ep heating, or to residual low-level ac- 
cretio n in the low state dStockman et al.|[T994l; IGansicke et al.l 
[T995h . In some polars, the geometry is favorable and the 
heated pole cap is eclipsed by the body of the white dwarf 
for part of the orbital cycle, allowing an accu rate determina- 
tion o f Teff from phase-resolved spectroscopy dGansicke et al.l 
l2006h . If the pole cap contributes at all phases to the UV 
light, or only phase-averaged data is a vailable, the data can 
be fit with a two-component m odel (e.g. I Gansicke et aLll2000l ; 
I Araujo-Betancor et al.ll2005bl) . T^ff from such analyses pro- 
vides an upper limit to the true WD temperature. 

4.4. Reliable T^-ff measurements 

In the context of using Teff as a measurable quantity that 
allows insight into the secular averages of the mass transfer 
rates in CVs, we have included in Table[T] only those sys- 
tems which we feel have a reliable Teff determination. Con- 
sequently, we omitted systems with published WD temper- 
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atures where the evidence for seeing the WD is ambiguous. 
Examples of such cases are WX Hyi, SSCyg, and RUPeg 
jSion & UrbanI 120021; iLong et al.l l2005h . where a plausible 
WD model fit to the UV spectra can be achieved, but no clear 
WD features are discerned (broad Lyman lines, narrow metal 
absorption lines). The decision of admitting a system to Ta- 
ble[T]is necessarily subject to a gray area, where some spec- 
troscopic evidence for the WD is present, but not sufficient for 
an accurate T,ff dete r mination, such as e.g . the ca se of Z Cam 
dHartlev et allfeoOSl) . iGan sicke & Koested d 19991) have shown 
that in the case of low spectral resolution and low signal-to- 
noise the UV data may be equally well described by a mod- 
erately hot WD or by an optically thick accretion disk. In 
this particular case, AH Men, the WD case can be excluded 
on the basis of the optical properties of the star, providing a 
clear warning against interpreting a slight flux turnover below 
1300 A as broad La from a WD photosphere. 

In the case of eclipse light curve analyses, we excluded 
a number of systems where we considered the data of too 
low quality, i.e. an at best marginal detection of the WD 
ingress/egress, as well as studies using oversimplified mod- 
els, such as approximating the WD emission in the different 
observed wave bands by blackbody radiation. 

In recent years, a number of strongly magnetic close 
WDh-MS binaries were identified that contain very cool 
(Teff S 9000 K) WDs and have mass transfer rates, as deter- 
mined from X-ray obse r vations, of a few 10~'^Myr~' (e.g. 
Reimers & Hagen '2000; ' Szkody et al.l l2003al: ISchmidt et all 
2005; Vogel et al. 2007). Given the fact tiiat many of them 
have MS companions that have spectral types too late to 
be Roche-lobe filling at the orbita l periods of the binaries, 
IWebbink & Wickramasingh^ (120051) suggested that these sys- 
tems are pre-CVs that have not yet evolved into a semi- 
detached configuration. The low mass transfer rates are com- 
patible with wind accretion, and the low WD temperatures 
mat ch with the predictions for th e average life time of pre- 
CVs lSchreiber & Gansickel(l2003h . Consequently, we exclude 
those systems from the present discussion. 

A final note concerns the intermediate polars (IPs), a class 
of weakly magnetic CVs in which the WD spin period is 
shorter than the orbital period, and partial accretion disks 
may form. In most IPs, the accretion rate is t oo high to dis- 
cern the WD even at UV w avelengths (e.g. iMouchet et al.l 
ll991l:lBeuermann et al.l2004l) . In a few IPs, moderately broad 
Balmer abs orption lines were detected in the optical (e.g. 
iHaberl et~ai.. 2002) but a more detailed analysis of the system 
parameters, in particular the dis tance, ruled out a WD photo- 
spheric origin of these features (de Martino et al."2006) In the 
case of EXHya (Eisenbart et al. 2002; Befle et al. 2003) and 
AEAqr (Eracleous et al. 1994) HST spectroscopy provided 
more convincing evidence for the detection of thermal emis- 
sion from the WD photosphere, compatible with temperatures 
~ 25 000 K, however, at least in the case of AE Aqr that tem- 
perature was cleary not that of the quiescent WD, but of the 
accretion-heated pole-cap. In summary, we did not include 
any IP in Table[T] because of the lack of a clear detection of 
the quiescent WD in any of these objects. 

5. DISCUSSION 

While we will now inspect the values listed in Table[T]for 
possible correlations, we must be utterly aware of the fact 
that the set of known CVWD temperatures is subject to se- 
vere selection effects. A very obvious, but crucial statement 
is that we need to be able to see the WD in order to measure 



its temperature, and, as mentioned above, this is not, or only 
marginally the case in systems where the mass transfer rate is 
too high. Therefore, it is possible that the temperatures ob- 
tained are rather lower limits than average values, as WDs in 
systems with higher accretion rates will be hotter - but not 
visible. This is most likely a stronger bias above the orbital 
period gap than below. 

Our collection of secure Teff measurements as listed in Table 
1 are plotted against system Porb in Figure |5] Non-magnetic 
(solid circles) and strongly magnetic (open diamonds) sys- 
tems are differentiated by symbol type. For the highest qual- 
ity measurements, the uncertainty in T^ff, which is dominated 
by the unknown WD mass, is indicated with an error bar 
Note that this means points without error bars are less accu- 
rate than those with, but their uncertainty is not easily quanti- 
fied. An approximate (M) sca le is shown on the right if an ex- 
pecte d typical mass (O.75M0: [Smith & Dhillonlll998l:lKnigg3 
|2006|) is chosen, but the degree of scatter expected is essen- 
tially u nknown. We sho w for comparison the empirical rela- 
tion of Patt ersoiil (Il984l) . using an assumed primary mass of 
M = Q.ISMq. From our Teff data we infer modestly higher 
mass transfer rates at all orbital periods, although a higher av- 
erage mass is also a viable explanation. 

It is useful to dislay a variety of theoretical predictions to 
compare with our data. Solid lines show the Tgff range ex- 
pected for the (M) obtained from typical interrupted mag- 
netic br^ing_scenarios, for the mass range 0.6 to O.9M0 
(see iHoweU et al.l 1200 ll and references therein). Here we 
have use d the mass-radius rela tion for the donor star in the 
study bv lKolb & Baraff j (1199^) for below the period gap and 
above the gap a value of (M) = 1O~''M0 at Porb = 3 hours and 
10~^Mo yr"' at forb = 6 hours, with a line between. Also 
shown is a cu r ve (do t-dashed) for the M = QJMp, history frorn 
iHowefl et all (|2001l) . The relation of lAndronov et alJ (l2003h 
for an unevolved donor is also shown, but implies much lower 
(M) than the data support. The mass-transfe r history of a case 
utilizin g the weaker braking law posited by llvanova & Taanj 
(l2004l) is shown by the dotted line. This curve corresponds 
to their (M) history for a 0.8Mq donor and a Q.6Mq primary, 
but for display purposes we have used M = Q.ISMq for the 
primary mass to obtain the Teff, thus the true predicted T^ff 
is likekly slightly higher For comparison to the magnetic 
systems, we have used the main sequence mass-radius rela- 
tion given by lHowell et al.l (1200 L) to determine (M) for non- 
conservative mass transfer under gravitational radiation only. 
The resulting T^ff is shown by dashed lines, again for the mass 
range 0.6 to 0.9Mo. 

Uncertainty and caveats due to selection bias will be dis- 
cussed below, as will the contrast between magnetic and non- 
magnetic systems. Comparing our measurements for non- 
magnetic systems to the predictions shows that we infer mass 
transfer rates that are lower than the predictions of traditional 
magnetic braking. Assuming that there is no strong selection 
bias toward lower (M) systems, not a trivial assumption as 
discussed below, our data above the period gap are more con- 
sis tent with the w eakened angular momentum loss proposed 
bv llvanova & T aam (2004), though still a bit higher than their 
predictions. The angular momentum law in the traditional 
magnetic braking picture was largely calibrated in order to 
raise the radius of the companion enough to reproduce the 
period gap, and, by virtue of this fact, evolution following 
the Ivanova & Taam (2004) relation might lack an appropri- 
ate period gap. The consistency of our data with the "softer" 
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TABLE 1 

Reliable T^ff measurements for WDs in CVs. Also listed are the CV subtype, the orbital period, the distance implied by the WD 

FIT, THE WD mass IF MEASURED INDEPENDENTLY FROM THE SPECTRAL FIT, AND THE DISTANCE OBTAINED FROM A TRIGONOMETRIC PARALLAX. 



System 


Type 


PcA [h] 


reff[K] 


± 


d[pc] 


M„d [Mq] 


d[pc] 


Ref 


GWLib 


DNAVZ 


1.280 


14700 




150-170 




,f^A 1 30 

104 ±30 


1,2 


BWScl 


DN? 


1.304 


14800 


900 


131±18 




3 


LLAnd 


DNAVZ 


1.321 


14300 


1000 


760 ± 100 






4 


EFEri 


AM 


1.350 


9500 


500 


~ 130 




163 ±«^ 


5,6, 2 


SDSSJ1610-0102 


DN? 


1.34 


14500 


1500 






7 


HS233 1+3905 


DN 


1.351 


11500 


750 


95 ±15 






8 


ALCom 


DN/WZ 


1.361 


16300 


1000 


800 ± 150 






9 


WZ Sge 


UN/ WZ 


1.361 


14900 


250 




U.OJ ih U.U4 


44 ± z 


lU, 1 1, IZ, 13,14, I 


^ W UiVid 




1 If^A 
1 . JD^ 


1 "^Qnn 

LjyvJyj 


yvAJ 








J 




DN'^ 


1 

i . JOO 


10500 


1000 








15 16 


HV Vir 


DNAVZ 


1.370 


13300 


800 


480 lb 70 




460 ± "° 


17 2 


WXCet 


DNAVZ 


1.399 


13500 




133 




18 


EG Cnc 


DNAVZ 


1.410 


12300 


700 


420 ± 60 






15, 2 


XZ Eri 


DN/SU 


1.468 


15000 


1500 




0.767 ±0.018 




19 


DP Leo 


AM 


1.497 


13500 




400 


0.6 




20 


V347 Pav 


AM 


1.501 


11800 


600 


1 77 + 33 

^ ' ' ^ 38 






21 


Rr" T TMn 


DN/SU 


1 503 


15200 


1 000 


285 lb 42 








EKTrA 


DN/SU 


1.509 


18000 


1200 


200 






22 


VY Aqr 


DN/WZ 


1.514 


14500 




187 




97 + " 


18 2 


OY Car 


DN/SU 


1515 


15000 


2000 






23 24 25 


VV Pnn 
V V r up 






1 1 Qnn 


ouu 


1 s 1 -\- 
151 ± 

144 ± 23 






Zi 


V534cen 


AM 


1.692 


14300 


900 






zl 


HTCas 


DN/SU 


1.768 


14000 


1000 






26,27,28 


VWHyi 


DN/SU 


1.783 


20000 


1000 




71 + 

^ 0.26 




29,30,31 


CUVel 


DN/SU 


1.88 


18500 


1500 






32 


MRSer 


AM 


1.891 


14200 


900 


, 1 S 

160 ±26 






21 


BLHyi 


AM 


1.894 


13300 


900 


163 






21 


ST LMi 


AM 


i.oyo 


10800 




1 1 S -k 99 
i i J ^ zz 






21 


EE Peg 


DN/WZ 


2.00: 


16600 


1000 


380 ± 60 






4 


DVUMa 


DN/SU 


2.138 


20000 


1500 




1.041 ±0.024 




19 


HUAqr 


AM 


2.084 


14000 










33 


QSTel 


AM 


2.332 


17500 


1500 








34 


SDSSJ1702+3229 


DN/SU 


2.402 


17000 


500 


440 ± 30 


0.94 ±0.01 




35 


AM Her 


AM 


3.094 


19800 


700 


90 






36,37,2 


MVLyr 


NL/VY 


3.176 


47000 




505 ± 50 






38 


DWUMa 


NL/VY 


3.279 


50000 


1000 


590 ± 100 


0.77 ±0.07 




39 


TTAri 


NL/VY 


3.301 


39000 




335 ±50 






40 


VI 043 Cen 


AM 


4.190 


15000 




200 






41 


WW Get 


DN 


4.220 


26000 


1000 








42 


UGem 


DNAJG 


4.246 


30000 


1000 




~ 1.1 


104 ±4 


43,44,45,13 


SSAur 


DN/UG 


4.391 


27000 










46 


V895 Cen 


AM 


4.765 


14000 


900 


511±- 






21 


RXAnd 


DN/ZC 


5.037 


34000 


1000 








47 



(l):Szkodv et al. 2002a, (2) Thorstensen 2003, (3) Gansicke et al. 2005, (4) Howell et al. 2002, (5) Szkody et al. 2006, (6) Beuennann et al. 2000, (7) Szkody et aU 
|5o67, (8) Arauio-Betancor et al. 2005a, (9) Szkody et al. 2003b, (10) Sion et al. 1995a, (1 1) Steeghs et al. 2001, ( 1 2) Long et al. 2004, (13) Harrison et aU 
12004, (14)SteeghsetaL 2007, (15) Southworth et al. 2006, (16) Littlefair et al. 2006b, (17) Szkody et al. 2002b, (18) Sion et al. 2003, (19) Feline et al. 2004 
(20) Schwope et al. 2002, (21) Arauio-Betancor et al. 2005b, (22) Gansicke et al. 2001, (23) Hessman et al. 1989. (24) Home et al. 1994, (25) Cheng et aU 
|2000, (26) Wood et al. 1992, (27) Wood et al. 1995, (28) Feline et al. 2005, (29) Gansicke & Beuermann 1996b, (30) Sion et al. 1996, (31) Smith et al. 200a 
(32) Gansicke & Koester 1999, (33)Gansicke 1999, (34) Rosen et al. 2001. (35) Littlefair et al. 2006a. (36) Gansicke et al. 1995, (37) Gansicke et al. 200a 
(38) Hoard et al. 2004, (39) Arauio-Betancor et al. 2003, (40) Gansicke et al. 1999, (4 1 ) iGansicke et al J200a (42) iGodon et aH2006L (43) lLong & Gilinandir999l 
(44) .Long et al..2006, (45) .Sion et aL1998, (46) .Sion et al..2004, (47) .Sion et al..2001i 
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Orbital Period (hours) 

Fig. 5. — Observed T^ff in the best observational cases (data from Table [T). Non-magnetic systems, a mix of dwarf novae and novalikes in low states, are 
indicated with solid circles and magnetic systems, all Polars, with open diamonds. Error bars indicate uncertainty dominated by the unknown WD mass. Points 
without error bars are less accurate than those with, but their uncertainty is not easily quantified. An approxmate mapping to (M) is shown on the right vertical 
scale assuming M = 0.75Mq , 0.6Mq or .9Mp, as indicated . The empirical relation of Patterson (1984; thick grey hne) is shown, along with several theoretical 
predictions: traditional magnetic braking (Howell et al.ll2001t between solid lines, dot-dashed line), Andronov et al. (2003; dot-dot-dash line), Ivanova & Taam 
(2004; dotted line), and evolution under only gravitational radiation losses (between dashed lines). There is a clear demonstration that, for Porb > 3 hours, stellar 
wind angular momentum loss is inhibited in the strongly magnetic systems. 
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law of llvanova & TaarnI ( |2004|) and the presence of the high- 
(M) VY Scl systems, suggest a very different picture, where 
the bloating of the companion might be localized to just above 
the period gap. There remains much room for improvement 
in the development of braking laws, as current laws remain 
fairly empirical with only modest input from understanding 
of the magnetospheric structure and even less from possible 
properties of the stellar dynamo. 

5.1. Uncertainty, Scatter, and Bias 

We would like to utilize our collection of Teff measurements 
mainly to constrain the dependence of (M) on Porb and thereby 
evaluate angular momentum loss laws utilized to predict this 
relation in CV evolution. There are two important sources of 
general uncertainty bearing on conclusions drawn from this 
dataset, which follow from the discussions of previous sec- 
tions: (1) uncertainty due to the effects of long-term M varia- 
tions on Teff, and (2) uncertainty due to inaccessibility of qui- 
escent Teff measurements at a given Porb- Each of these will 
introduce qualifications to the naive interpretation that Lq can 

be converted directly into (M) representative of that Poib in- 
terval for that kind of system (magnetic or non-magnetic) via 
a relation like equation [l] 

5.1.1. Long-term M Variations 

Given only a single T^ff measurement, it is always possible 
that (M) is something other than that implied by T^ff, and we 
are observing a transient state. However, we do not have only 
one Teff measurement; we have multiple measurements for ob- 
jects in each of the 3 classes discussed in section lU dwarf 
novae, novalikes displaying low states, and polars. In light of 
the time series and analysis of section |2l medium-term varia- 
tions in M, on scales of hundreds to thousands of years, would 
manifest as object-to-object scatter in T^ff at a given Porb- In 
contrast, the measurements of T^ff are remarkably consistent 
for the groups. The novalikes show the best evidence of scat- 
ter, and therefore of M variations on these timescales. They 
also would have the shortest adjustment time for their Tgff due 
to the large implied (M) and luminosities. Thus novalikes 
displaying low states might accrete at ~ IQ~^Mq yr"' for pe- 
riods of 100 years or so, and but have (M) of more like a few 
lO-'^M© yr-i. This can be better quantified as more measure- 
ments become available. 

The consistency of the T^ff values at a given PqA within 
each group provides good evidence against large, medium- 
term variations in M. For the (M) appropriate for objects 
with Teff < 20kK, this statement extends to even 10^ years 
or more. That is, the objects for which we have good mea- 
surements of Teff appear to have remarkably similar M histo- 
ries over the last 10^ to 10^ years, though these objects were 
likely to have been born at a variety of Porb's. As noted ear- 
lier, a true scatter is expected due to variations in Mace and M 
among observed systems. This makes the consistency among 
the measurements even more remarkable; enough to suggest 
there might be some mechanism, possibly in Classical Nova 
outbursts, which regulates M beyond just selection bias. More 
careful treatment of sample bias, in particular with respect to 
the WD mass, which actually entails more uniform UV study 
of known systems, would be necessary to make conclusions 
on the presence or absence of such a mechanism. 

The extreme case of uncertainty due to long term M varia- 
tions arises from the fact that inactive systems, as e.g. in the 



hibernation scenario of lSharaet all (11981) or in irradiation- 
induced mass transfer cycles dRitter et al may not be 
included in the CV census at all, and therefore there is no op- 
portunity to measure their T^ff. This can be roughly quantified 
by considering a system with a duty cycle /, which is active 
for a time period factive and then ceases mass transfer for a pe- 
riod of ractive(l -/)//■ During the active phase, the Lq of such 
a system is given approximately by 

iq^acdve ~ Lq({M)) + (5Lq(Mactive , factive) (6) 

where the dependence of SLq on factive is the same as the de- 
pendence on thermal time shown in Figure |2] We introduce a 
fractional response function Rl such that (5Lq(Mactive, factive) = 

iq(Mactive) ' ^L(f active) SUCh that < Pl < 1- Then PL(factive) 

is the fractional response to an order unity variation in M 
on a timescale factive, which increases with factive, and is the 
unitless quantity actually shown in Figure |2] Noting that 
Mactive « {M) I f and that Lq oc M, we have 

iq.active ~ L^({M)) ■ (1 +-R£(faetive)//) • (7) 

This demonstrates that for short factive, such that P^ <C 1, 
i'q.active providcs a good proxy for Lq((M)), and therefore Leff 
is a good indicator of (M) even with a direct conversion. 
However, for PL(factive) ^ 1 and small /, Lq during the ac- 
tive state can be largely unrelated to (M), and instead be set 
by / ~ (M) /Mactive and factive- It should be noted that a long 
factive and a low / can imply a very long recurrence time. 

The low duty cycle and long recurrence time scenario just 
described is very unlikely to apply to systems with Leff < 
20kK. The proximity of the (M) implied for these systems to 
the lower limit set by gravitational radiation angular momen- 
tum losses (dashed lines in figure |5]l excludes / ^ 1 . This 
means that time series like those presented in Section 12.41 
where / = 0.5, are appropriate for these systems. Addition- 
ally, as discussed above, very long factive <^ 10^ years would 
be necessary to ensure consistency among so many indepen- 
dent objects. 

The situation for non-magnetic systems above the period 
gap is less constraining. Assuming that there is some mech- 
anism which can regulate Mactive with some precision, the 
consistency of several measurements in this region implies 
that factive must at least be a few times the time it takes to 
reach Lq(Mactive)- From the higher (M) curve in Figure |2] it 
will take approximately 5,000 years to rise within 20% of the 
Lq (Mactive) indicated for these systems, so that we can esti- 
mate factive ^ 10"^ ycars. In this case we are assuming that 
P/,(factive) — 1, so that in order to overestimate the (M) by a 
factor of 10 requires / = 1/10, and thus a recurrence time of 
10^ years and a similar duration of the inactive phase. Thus 
long-timescale hibernation scenarios (e.g. S hara et al.lll986l) 
cannot be excluded by the current Leff data. There is also 
no apparent evidence favoring such scenarios, in the form of 
downward scatter of objects still transiting between inactive 
and active phases. 

A sample of Left measurements for detached systems in the 
3-6 hour Poib interval could conclusively rule this out by the 
absence of excess WDs with Leff ~ 15 kK compared to longer 
periods. Such samples are currently being constructed and 
the results are ambiguous. Several detached WDh-MS sys- 
tems that bear the characteristics of hibernating CVs have 
been identified, namely BPM 71214 with Ppfb = 290 min an d 
Leff = 17000K, (iKawka et alJl200l iKawka & Vennesll2003h . 
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EC 13471-1258 with f pi-b = 217 min and T^n = 14220 ± 350K 
(lO'Donoghue et alJ l2003h . and HS 2237+8154 w ith forb = 
178 min and Teff = 1 1 500 K jGansicke et al.ll2004a l). The first 
two show promise, however there is significant selection bias 
toward finding hot WDs and there is an expected (contami- 
nant) population of systems which are just coming into con- 
tact and simply have young WD primaries. 

5.1.2. Inaccessible Quiescent T^fi Values 

The second major source of uncertainty in drawing con- 
clusions from the available set of Teff measurements is due 
to the set of circumstances which must come to pass in or- 
der to allow direct measurement of the WD photosphere. For 
Dwarf Novae, which accrete in bursts with M (M), it is 
also necessary to wait a sufficient period of time after the out- 
burst in order to get a good idea of the baseline quiescent T^ff 
that is escaping from the deeper layers of the envelope with 
longer thermal times. As mentioned above, this latter can be 
achieved largely empirically by selecting the timing of Teff 
measurements with respect to disk outbursts. 

Non-magnetic Cataclysmic var iables with Ppfb < 2 hours 
are predominantly Dwarf Novae (iRitter & Kolbll2003h . The 
selection criteria for such systems are important: the emission 
of unknown origin discussed in section l^Tl must be less bright 
than the WD in the UV, and the absorption of the system must 
be low enough that the WD can be measured well, placing 
a constraint on the inclination of the system. Without better 
characterization of the unidentified broad-band emission, our 
only option is to assume it is a random contaminant which is 
uncorrected with (M). Again, it is possible to confirm this 
with better UV study of known systems. Since both of these 
selection criteria are not expected to correlate with (M), we 
believe that our sample of non-magnetic systems in the forb < 
2 hour range should be representative of the typical (M) in 
these systems. There might be a slight bias toward low (M) 
due to their having more accessible quiescent intervals, but 
there are no indications that this is the case. 

For forb > 3 hours, DN are a minority of the population, 
and thus there is concern that only particular systems with low 
(M) have made it into our sample. If true, this would imply 
that the region at higher T^ff than the measured systems in the 
3.5-5 hour period range in Figure |5] should contain the higher 
(M) systems which did not enter our sample. This would im- 
ply a higher (M) than we can infer directly from the measured 
systems, and thus impart more favor to the traditional mag- 
netic braking prescriptions. 

5.2. VY Sculptoris stars 

The three VY Scl stars with well-determined temperatures 
stand out in the 3^h period range containing the hottest 
CVWDs known, and consequently have very high mass trans- 
fer rates. In effect, the deduced mass transfer rates exceed 
those predicted by the standard evol ution theory for the major- 
ity of CVs within that period range (lKolblll993tlHowell et al.1 
1200 lb . Teff measurements of VY Scl require the fairly prompt 
observational attention once they enter a low state, preferably 
with an UV facility, which explains the small number of avail- 
able values. Practically all VYS cl stars are located withi n 
the 3-4 h orbital per iod range (e.g.THoneycutt & Kafkal2004l) . 
Furthermore, Rodrig uez-Gil et aTl (l2007h have shown that the 
SWSex stars, intrinsically bright novalike variables which 
are intimately relate d to the VY Scl stars (.Hameury & Lasotai 
l2002t lHeUieiil2bO0t in fact, the two groups overlap to a large 



extent) are the dominant population of CVs in the 3^ h or- 
bital period range. Speculating that high T^ff and (M) are a 
common characteristic to all VYScl/SWSex stars suggests 
that these systems represent an exceptional phase in CV evo- 
lution. One possible explanation is that these are systems that 
just evolved into a semi-detached configuration, as the the 
mass transfer goes thr ough a short peak during turn-on (e.g. 
iD'Antona et al.llT989f) . and that CVs are preferentially born 
within the 3-4 h period range, which would be the case if the 
initial mass distribution is peaked to wards equal rti asses in the 
progenitor main-sequence binaries dde Koollll992h . 

5.3. Polars and Wind Braking 

In agreement with the interrupted magnetic braking sce- 
nario for non-magnetic CV evolution, these T^ff measurements 
indicate that (M) is approximately an order of magnitude 
larger above the period gap than below. But Teff measurements 
can do better than this. Knowing the relation between T^ff and 
(M) we can say that the objects below the gap are roughly 
consistent with gravitational radiation losses, with possibly 
some enhancement of a factor of 2 or 3 for non-magnetic sys- 
tems. In contrast (M) above the period gap is an order of 
magnitude greater than that predicted by gravitational radia- 
tion alone. There is, finally, an additional constraint that is en- 
tirely specific to the magnetic braking mechanism: we find a 
marked difference between the (M) implied for non-magnetic 
systems and those of strongly magnetic systems (Polars) and, 
as shown in Figure|5] the (M) in Polars is consistent with that 
expected from gravitational radiation angular momentum loss 
alone, while that in non-magnetic systems is at least an order 
of magnitude higher As shown in section 13.21 such a decre- 
ment in Teff, which is measured away from the poles where 
the accretion impacts, is too large to be explained by accre- 
tion geometry for the magnetic fields observed. Additionally, 
from the discussion in section 15.1.11 such a large difference 
would require an extreme assumption about the duty cycle in 
non-magnetic systems. Therefore the contrast between mag- 
netic and non-magnetic systems must arise from a difference 
in (M). 

The lack of an enhanced (M) in magnetic systems arises 
from changes in the magnetic field structure near the sec- 
ondar y which hinders the loss of angular momentum via a 
wind dLi et al.| [r994b). Our measurements provide the best 
direct evidence that this does occur, and additionally that the 
resulting (M) is consistent with gravitational radiation. This 
provides very strong support for the basic picture of magnetic 
braking, though the precise mechanism by which it ceases is 
still somewhat mysterious. We should highlight that this re- 
duction of magnetic braking in Polars is widely expected, and 
was originally proposed to e xplain the lack of an ap parent 
period gap in rnagnetic CVs dWick ramasinghe & Wu 1994t 
iLi et al.lll994i IWebbink & Wickram asinghe 2002) . The re- 
sulting slower evolution of magnetic systems will also en- 
hance the number of magnetic CVs relative to nonmagnetic 
ones with respect to field WDs, givi ng roughly the fraction 
observed (iTownsley & Bildstenll2^5h . 

5.4. Below the Period Gap 

Due to the need for low (M), many of of the high-quality 
Teff measurements are for systems with PqA < 2 hours. An 
expanded plot of this region is shown in Figure |6] For M = 
0.6Mq we use the (M) from lKolb & Baraffd (Il999h directly 
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Fig. 6. — Closeup of Pojb < 2 hours, filled circles indicate non-magnetic 
systems and open diamonds indicate magnetic systems (see Figure[5j- Over- 
plotted are the range of T^ff expected for accretion driven by gravitational 
radiation tKolb & Baratfe 1999 ) for two WD masses. Note that the theoreti- 
cal minimum period is notably shorter than that inferred from the data. 

and for I.OMq we use their M-R relation for the secondary 
and gravitational radiation angular momentum losses. We 
display the predicted region for O.OSMign < M^^c < 0.95Mig„ 
which object will traverse during buildup toward a classical 
nova when Mace = Mgn- Many of the measurements for non- 
magnetic objects are clustered about Tetf =15 kK, providing 
evidence, as discussed above, against long-term M variations. 
The exception to this is that there is significant downward 
scatter at the shortest orbital periods 1.3-1.4 hours. The pres- 
ence of this scatter at these forb and not at slightly longer 
ones strongly suggests that the (M) in these objects has begun 
to decl ine as expected wh en they pass beyond period min- 
imum (iHowell et al.ll200Tl) . This feature highlights the dis- 
agreement of the observed period minimum and that predicted 
by theory (Kolb & Baraffe 1999), which is apparent in Fig- 
ure |6] because the minimum is not very sensitive to M. As 



discussed in iTownslev & BildstenI (l2003l) we find that either 
short-period CVs have M Q.9-l . QMp^ (s uch as e.g. found in 
SDSS 1035-H0555. [Littlefair et all (l2006bl) ) or mass transfer is 
enhanced by a factor of 2-3 over that predicted by only gravi- 
tational radiation losses. There is also some evidence from the 
few measurements near Porb = 1.8 hours that the dependence 
of (M) on Porb is not as flat as is predicted by gravitational 
radiation. More measurements will be necessary before con- 
clusions can be drawn. 

Even here below the period gap there is a smaller but signif- 
icant difference between strongly magnetic and non-magnetic 
systems. The energy released by compression deeper than 
the point at which the material is able to spread over the star 
amounts to better than 80% of the Lq expected in the non- 
magnetic case. However such a decrease in Lq only reduces 
Teff from 15 kK to 14 kK, not enough to account for the dif- 
ference in Tgff between magnetic and non-magnetic systems 
below the gap. For high-field cases, Z? ~ 10** G, the reduction 
can reach 50%, bringing Teff down to 1 3 kK. However, none of 
the systems used in this paper (Table[T]i has such a high field 
(see Sect. l4.3l ). It still appears that either (M) is slightly above 
that due to gravitational radiation or these magnetic objects 
do have typical M ~ O.IMq. A typical mass above 0.6Mq 
is expected both from simple selection bias because the lumi- 
nosity increases with M and since field magnetic WDs tend to 
be higher mass than non-magnetics. 
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APPENDIX 

SIMULATION OF HYDROSTATIC ENVELOPE IN PLANE PARALLEL 

In this appendix we detail our treatment via simulation of the envelope of the accreting WD. We work in the plane parallel, 
hydrostatic approximation so that P = gy at every point. Discretizing the temperature in y and time, the heat equation, Q, is 
integrated forward in time using a Crank-Nicholson like integration rule 
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' dT 
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dt 
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where the derivatives of T in y on the grid are evaluated via centered differences 
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Initial conditions are taken from the static models used by TB04. Note that eq. ( lAlb is an implicit integration rule, representing 
equations in unknowns, which is solved for T^""^' using a Newton-Raphson iteration. Timesteps are chosen such that the 
maximum AT /T at any point on the grid is less that 10"^. The steady state solution produced by the time-dependent code 
compares very well with the static solutions produced by direct integration of the structure equations in TB04. 

The flux at the outer boundary is found by integrating eq. (O in from the photosphere using a 4th order Runge-Kutta integrator 
A shooting method (root-find) is used to vary /^photosphere to match Tiyo^er) at the outer edge of the simulation grid, jouter^ and the 
resulting /^(jouter) is used for the flux at the edge of the grid. While there are no convection zones in the simulation domain, there 
is a convec tion zone in the part of the envelope which forms the boundary condition. Convection here is treated with the ML2 
formalism (iBergeron et al.|[l992h . Use of this boundary condition is equivalent to the approximation that the thermal time of the 
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layer of depth jouter is negligible. The inner boundary condition is constant temperature, fixed to the equilibrium temperature at 
the (M) under consideration (TB04). 

Two different grids are used in this work. For the long-timescale variation studies of section l24l we use = 120 points evenly 
spaced in logy extending from fromy = 6 x 10'' g cm"^ to 3.9 x 10'^ g cm~^, corresponding to mass coordinates of 1.6 x lO'^M© 
and O.OIMq, respectively on a 0.9Mq WD, R = 6.45 x 10** cm. For (M) = 5 x 1O""M0 yr"' at this mass, the thermal timescale 
of the outermost point fth(3'outei) — 10 days, much shorter than any variability being studied. For the study of TT Ari in section[3] 
the upper edge of the grid was extended to shorter thermal timescales. In this case we used 100 points between y = 2.75 x 10^ 
and 3.64 x 10** g cm"^ and 50 points extending down to 3.9 x 10'^ g cm"^, for a total 150. This outer point corresponds to 
fthCyouter) — 15 hours for the parameters used to reproduce TT Ari. 

The depth of the boundary between the solar abundance accreted material and the underlying 50/50 carbon and oxygen material, 
jHHe = A/acc/(47r/?^) is tracked with a separate variable whose evolution is directly specified from M{t). Although selecting 
between two compositions depending upon whether y^ < yuHe or ynne < yk gives suitably accurate results, a smooth progression 
is significantly more easily integrated numerically. This allows a smooth change in abundance at a given point while still allowing 
a sharp interface to be represented on a modest resolution grid. We now proceed to describe our interface treatment in detail. 

Consider the interface as being at column depth yi = ymin lying between two gridpoints at depths yk and yk+\ ■ The temperature 
at each of these points is 7i and Tk+\ respectively, and that at the interface is T/. Both the temperature and flux must be continuous 
at the interface, so letting Fh and Fc indicate the flux evaluated at yi on the H/He side and the C/O side respectively, we have 



AacT^ dT 

3k dy 
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Kc(yi,Ti) \ dy 
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where the subscripts H and C on the derivatives indicate evaluation on respective sides of the interface and on k indicate evaluation 
at Pi and 7} but with H/He and C/O composition respectively. It is inadvisable to attempt to solve for Tj directly because yi-yk 
can become arbitrarily small and lead to singularities. Instead we wiU solve for the derivatives. The interface temperature can be 
written by expanding from both directions 



7} = Tk + (y,-yk) 
which upon combination with eq. ( IA3b gives 
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kh(Ti) 

By using the first expression for Tj from ( IA4l l. this can be solved for (dT/dy)H- 

Finally, all the derivatives near the interface can then be constructed from the grid quantities and (dT /dy)H in a way which 
mimics centered differencing. We let (^)h stand in for the first order difference at the midpoint between k and k+ 1, since it 
should be approximately what that first order difference would have been if there were no change in composition. This gives 
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For the outer portions of the envelope we use the 2002 update to the OPAL equation of state tables dRogers et alj|1996h . and for 
higher densities we use the analytical approximations for a fully ionized plasma from lPaczvriskil (119831) and Coulomb correction 
from Chabrier & PotekhinI ([199 81) . While these two EOS methods are very consistent at the table edge, a linear average in a 
crossover region of a factor of 5 in density i s used to smooth the boundary. OPAL radiative opacities (.Iglesias & Rogersll 19961) 
are also used along with conductivities from lltoh et alJ (Il983l) . 
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